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Goal  programing  (6?)  or,  more  appropr lately  generalized  goal  pro¬ 
graming  (GGP)  [21,  24,  27],  Is  on*  of  several  approaches  that  have  bean 
proposed  for  the  modeling,  solution  and  analysis  of  the  multiobjective 
mathematical  programming  problem.  However,  GP  has  distinguished  itself 
from  most  alternate  multiobjective  methods  in  its  computational  efficiency 
and  evidence  of  widespread,  actual  implementation.  Even  the  most  adamant 
of  its  critics  recognizes  that  GP  la,  as  of  now,  the  "workhorse"  of  the 
multiobjective  mathematical  programing  methods. 

The  computational  efficiency  of  GP  is  particularly  evident  in  dealing 
with  either  linear  or  nonlinear  multiobj active  models  having  continuous 
variables  [21,  24,  27,  28;  29].  Here,  GP  has  shown  equivalent  performance 
to  that  available  in  conventional  (i.e. ,  single  objective)  algorithms. 
However,  when  deal J  .g  with  models  which  involve  integer-valued  variables, 
the  performance  of  various  (exact)  integer  GP  algorithms  has  been  less  than 
impressive.  As  a  result,  the  majority  of  truly  large-scale  integer  GP 
problems  have  been  solved  by  various  heuristic  approaches  [20,  21,  23,  30]. 

Relatively  recently,  there  have  been  some  substantial  accomplishments 
in  solving  single-objective  integer  programming  problems  via  generalized 
networks.  For  some  problems,  such  an  approach  has  been  shown  to  be  up  to 
several  hundred  times  faster  (and  with  less  round-off  error  and  reduced 
storage  requirements)  than  available  through  conventional  software  [6,  7,  8, 
12,  13,  14,  19,  37,  39].  One  deficiency  of  such  approaches,  at  least  in 
the  opinion  of  the  multiobjective  advocate,  has  been  the  focus  on  a  single 
objective. 


I 

! 


2 

Purpose 

The  purpose  of  this  pspsr  is  to  rsport  on  sn  ongoing  rssnsreh  effort 
that  is  directed  toward  the  development  of  a  hybrid  approach  which  combines 
GP  and  Generalised  Networks  (GN) .  As  a  result,  we  term  the  hybrid  approach 
"GP-GN".  A  few  previous  papers  [22,  26,  36,  38]  have  addressed  earlier 
results  in  this  area  and  thus,  herein,  we  expand  on  this  topic  with  our 
emphasis  focused  primarily  on  the  modeling  aspects  of  GP-GN  since  there  is 
an  abundance  of  literature  readily  available  in  regard  to  the  algorithms  of 
generalized  networks  [1,  2,  3,  4,  5,  6,  7,  8,  9,  10,  11,  12,  13,  14,  IS,  16, 
17,  18,  19,  31,  35,  37,  39,  40]. 
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2.  TWO  CP  MODELS 


Thar*  ar*  a  wide  variety  of  GP  models  including  vaightad  GP,  lexico¬ 
graphic  GP  (i.a. ,  tha  uaa  of  the  so-called  "preemptive  priority"  concept), 
minaax  GP  (which  includes  fussy  prograaerimg) ,  and  Interactive  GP  (which  is 
used  to  generate  a  subset  of  nondoainated  solutions)  [24,  25,  27,  30,  34]. 
However,  in  this  paper  ve  shall  initially  restrict  our  attention  to  two 
particular  forms  of  GP:  (1)  weighted  GP  and  (2)  lexicographic  GP.  The 
mathematical  models  used  to  represent  these  two  types  of  GP  are  then 
presented  below. 


Weighted  GP  Model 

Here,  we  wish  to  find  x  so  as  to: 

y 

minimize  z  (u^  +  w^) 


(2:i> 


where: 


ft  (x) 


I  < 


\*l 


ieF 


f±  (x)  +  ieG 


x,  n,  P  >  0 


(2.2) 


(2.3) 

(2.4) 


F  •  the  set  of  rigid  constraints 

G  -  the  set  of  goals  (where  the  right-hand  side  value  represents  a 

target  value,  or  aspiration  level,  for  the  specific  goal) 

-  the  negative  deviation  of  goal  i 
p^  «  the  positive  deviation  of  goal  1 
ui>wi  -  the  weights  associated  with  and  respectively 
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Lexicographic  GP  Modal 

If  the  problem  goals  are  rank  ordered  according  to  importance,  we  can 
model  the  problem  as  follows.  Find  ic  so  aa  to  determine  the  lexicographic 
minimum  of: 

a  ■  (a^,  •••,  •  ••»  (2.5) 

s.t. 

fA  (x)  -  P±  -  b±  Vi  (2.6) 

x,  n,  P  >,  0  (2.7) 

where: 

a  *  the  "achievement  vector"  for  which  we  seek  the  lexicographic 
minimal  value. 

a^  -  a  function  of  the  goal  deviation  variables  (n,p)  that  are  to  be 
minimized  at  priority  level  k.  That  is: 

»k  *  gk  (n.p)  (2.8) 

Solving  the  Models 

If  the  decision  variables  (x)  are  permitted  to  take  on  fractional  values, 
computationally  efficient  algorithms  exist  that  may  be  used  to  develop  the 
solution.  For  example,  if  the  models  involve  strictly  linear  functions 
(i.a.,  all  f^  (x)  are  linear)  then  algorithms  have  been  developed  for  use 
on  IBM,  UNIVAC  and  CDC  computers  (as  veil  as  others)  that  solve  problems 
of  the  same  size  (and  with  approximately  the  same  efficiency)  as  available 
through  the  use  of  conventional  simplex  codes.  For  example,  the  solution 
of  lexicographic  linear  GP  models  with  up  to  16,000  rows  (i.e.,  goals  and 
rigid  constraints)  by  as  many  variables  as  may  be  stored  is  available  on 
IBM  360  or  370  computers  having  the  MPSX  software  package  [24,  27,  29]. 


Ia  the  um  of  nonlinear  (weighted  or  lexicographic)  GP  nod  ala  (again, 
with  continuoua  variables),  computer  codas  have  bean  developed  (and  applied 
Co  a  variety  of  real  world  applications  [21,  24,  28])  that  will  solve 
problans  with  several  hundred  rows  by  thousands  of  variables.  One  actual 
Implementation  (1978  [28])  Involved  the  solution  of  a  problem  with  400 
rows  and  10,000  variables  and  was  solved  in  approximately  nine  minutes  on 
the  CDC-7600  computer. 

However,  as  mentioned,  the  efficiency  of  computer  codes  for  GP  models 
involving  integer  variables  is  not  particularly  impressive.  To  the 
author's  knowledge,  meet  truly  large-scale  integer  GP  problems  actually 
dealt  with  have  been  solved  by  means  of  heuristic  techniques.  Further, 
such  heuristics  have  typically  been  carefully  tailored  to  handle  the 
specific  problem  under  consideration  and,  thus,  an  approach  to  the  general 


case  has  been  lacking. 


3.  INTEGER  PROGRAMMING  AND  NETFORM 


There  are  nuaerous  approaches  to  the  modeling  (and  solution)  of  gener¬ 
alized  networks.  For  purposes  of  illustration  only,  ve  will  employ  models 

and  notation  used  in  one  of  these  approaches,  known  as  NETFORM  [6,  7,  12,  14]. 
In  this  section,  we  shall  briefly  summarize  the  NETFORM  modeling  approach  to 
conventional  (i.e.,  single  objective)  problems. 

The  NETFORM  Arc 

The  building  blocks  of  the  NETFORM  model  are  its  nodes  -  and  the  arcs 
connecting  these  nodes.  Figure  3-1  depicts  a  general  NETFORM  arc  and  its 
associated  nodes.  The  origin  node  is  node  "i"  and  the  terminal  node  is 
node  "j".  The  (directed)  arc,  ”i-j",  connects  these  nodes.  The  "cost"  per 
unit  flow  is  given  as  c  .  and  is  shown  within  the  rectangle, 

J 


(L,U)* 


Figure  3-1:  NETFORM  Arc 


The  flow,  x.  . ,  has  a  lower  bound  (L)  and  upper  bound  (U)  as  indicated  in 
J 

the  parentheses.  Note  that  an  asterisk  will  appear  only  if  the  flow  across 

the  arc  must  be  integer  valued.  A  particularly  important,  and  powerful, 

feature  of  the  model  is  given  in  the  triangle  as  G  . .  G  .  is,  in  turn, 

i*j 

the  gain  factor  which  reflects  the  amplification  or  attenuation  of  the  flow 
through  the  arc. 
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One*  a  problem  has  baan  ao  modeled,  various  software  is  available  for 
its  solution  and,  as  mentioned,  the  results  can  be  extremely  Impressive. 

One  application,  for  the  U.  S.  Treasury,  Involved  a  model  with  5,000  nodes 
and  625,000  arcs.  The  problem  was  reportedly  solved  in  9  minutes  on  a 
UNIVAC  1108  at  a  cost  of  $90.  A  good  conmercial  LP  code  would  have,  under 
best  conditions,  required  about  20  hours  of  CPU  time,  at  a  cost  of  about 
$24,000  [14] . 

The  Conventional  Integer  Programming  Model 

The  generalized  network  or  NETF0RM  modeling  approach  is  not  restri' 
to  just  those  problems  which  naturally  take  on  a  network  representation 
(e.g.,  assignment,  transportation,  transshipment).  We  can,  for  example, 
transform  any  single  objective  (and.  In  the  next  section,  multiobjectlve) 
integer  programming  problem  into  an  equivalent  GN  model.  This  is  best  illus¬ 
trated  via  a  numerical  example.  Consider  the  following  model  in  which  all 
variables  must  cake  on  0  or  1  values: 

minimize  3x^  +  7x2  +  x^ 

S  a  £•  • 

(I)  xl  +  x2  +  x3  —  2 

(II)  -2xx  +  5x2  >3 

Xj  e  (0,1)  J-1,2,3 

The  equivalent  GN  model  will  contain  one  "objective"  node  and  this  will 
serve  as  the  source  node  for  the  network.  Also,  there  will  be  a  node,  for 
each  variable,  connected  directly  to  the  objective  (or  source)  node.  Finally, 
there  will  be  a  node  for  each  constraint  and  these  constraint  nodes  will  be 
directly  connected  to  those  variable  nodes  as  associated  with  the  variables 
which  appear  in  chat  constraint.  This  leads  to  the  GN  model  shown  in 
Figure  3-2. 


The  "availability"  of  tha  objactlva,  or  source,  ooda  la  shown  aa  <3 
that  up  to  thraa  varlablaa  may  taka  on  valuea  of  ona.  Tha  requira- 
aanta  at  tha  terminal,  or  constraint,  nodes  are  thoaa  associated  with  each 
constraint  (aaa  original  modal).  Tha  coat  of  variable  1  In  tha  objactlva 
is  3  units  as  shown  on  the  arc  from  tha  objective  node  to  the  variable  one 
(x^)  node.  Tha  "gain"  across  this  arc  is  2  units  which  means  that 
variable  one  appears  in  exactly  two  constraints.  The  arcs  from  the  variable 
nodes  to  their  associated  constraint  nodes  have  a  gain  factor  equal  to 
their  technological  coefficient  in  that  constraint.  Solving  the  GN 
model  of  Figure  3-2  will  then  provide  the  solution  to  the  original  mathe¬ 
matical  model. 
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4.  WEIGHTED  ®-GH 

We  shall  first  address  the  modeling  of  the  weighted  integer  GP  model 
via  GN.  The  resulting  concepts  and  approach  may  then  be  extended  so  as  to 
encompass  the  subject  of  the  following  section:  lexicographic  integer  GP. 

The  easiest  way  to  explain  the  GP-GN  approach  is  via  a  simple.  Illus¬ 
trative  example.  Consider  the  following  weighted  Integer  GP  model: 

■  Find  x  so  as  to 


minimize 

z  -  5n^  +  2p2 

(I) 

*1  + 

x2  +  1  2 

(II) 

3x! 

+  2x3  +  -  p1  -  4 

(III) 

3*1  + 

2x2  +  7x3  +  n2  “  p2  "  5 

Xj  6 

(0,1) 

Approach  One 

Recalling  the  modeling  approach  used  in  the  conventional  integer  pro 

grammlng  model,  as  Illustrated  in  the  previous  section,  we  note  that  it  may 
be  directly  applied  to  our  example  since  it  too  is  in  the  form  of  a  single 
objective  Integer  model.  This  approach  would  lead  to  a  GP-GN  model  as 
shown  in  Figure  4-1.  Here,  we  simply  treated  the  goal  deviation  variables 
(n^,  A 2%  P-j.  and  o^)  as  °>odel  variables  and  added  the  corresponding  variable 
nodes . 

We  choose,  however,  not  to  use  this  approach  but,  rather,  employ  an 
equivalent  alternate  modeling  scheme,  in  the  construction  of  the  weighted 
integer  GP-GN  model.  We  discuss  this  technique  below. 
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Approach  Two 

This  second  approach  involves  tha  incorporation  of  a  n aw,  "complex" 
nods.  This  special  nods,  denoted  as  tha  goal  nods,  is  coaposed  of  a  staple 
terminal  node  plus  two  additional  ares.  A  general  goal-node  is  shown  in 
Figure  4-2.  In  this  figure,  we  note  that  node  A  is,  in  essence,  a  "soft" 
constraint  node.  That  is,  it  allows  for  some  deviation  either  above  or 
below  the  goal  target  value,  b^. 

The  input  to  node  A  reflects  the  actual  value  achieved  for  that  goal, 
which  we  may  call  a^.  The  desired  value,  b^,  is  shown  to  the  right  and,  in 
the  network,  is  a  sink  node  requirement.  However,  if  a^  >  then  there  is 

positive  deviation  (or  "overflow")  which  flows  out  of  node  A  into  the  top, 
positive  deviation  (p^)  branch  at  a  cost  of  w^  units  par  unit  of  flow.  On 
the  other  hand,  if  ai  <  b^,  then  there  is  negative  deviation  (or  "underflow") 
which  requires  a  flow  into  node  A  -  as  depicted  by  the  lower  branch  (n1) . 

Tha  flow  on  this  branch  Incurs  a  cost  of  u^  units  per  unit  of  flow.  Such 
a  goal-node  will  then  be  associated  with  each  goal  In  the  GP  model  under 
consideration.  Combining  the  GN  modeling  process  with  the  goal-node  concept, 
we  may  model  our  example  as  shown  in  Figure  4-3.  (Recall  that  numbers  in 
rectangles  represent  coefficients  or,  in  this  case,  weights  while  these 
in  triangles  refer  to  tha  gain  across  the  arc.) 
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Figure  4-2:  The  Goal-Mode 


To  summarize,  Che  use  of  Che  goal-node  leads  Co  a  weighted  inceger  GP 
model  wherein: 

(1)  There  Is  a  single  source  node  (node  0)  wich  a  "supply"  equal  Co 
che  number  of  zero-one  variables  In  che  problem. 

(2)  There  are  arcs  from  che  source  node  directly  co  a  set  of  variable 
nodes  (i.e.,  decision  variables  only).  The  multiplier  (gain)  on 
Che  arc  indicates  Che  number  of  constraints  plus  goals  in  which 
che  associated  variables  appears. 

(3)  The  variable  nodes  are  connected,  via  arcs,  to  che  constraint  and 
goal-nodes  in  which  they  appear.  The  multiplier  for  these  arcs 
is  equal  to  the  technological  coefficient  of  the  variable  in  the 
respective  constraint  or  goal. 
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(4)  A  constraint  noda  (aaa  nods  I)  is  than  a  • lapis,  tarainatlng  nods. 
The  "demand"  of  that  noda  is  shown  to  its  right.  A  goal-node  is 
also  a  terminating  node  but  is,  as  noted,  a  complex  node  (as 
shown  for  goals  II  and  III) . 

When  the  standard  a s simp t ions  (see  the  references)  of  generalised  net¬ 
works  are  satisfied,  solution  to  this  form  of  ehe  weighted  Integer  GP-GN 
model  is  available  via  the  appropriate  generalised  network  software.  The 
solution  thus  found,  for  our  example,  would  be: 
x  *  -  (1,0,0)  or  (1,1,0) 

and  s*  •  5 

We  now  proceed  to  the  second  form  of  the  GP  model,  l.e.,  lexicographic  GP. 


16 


5.  LEXICOGRAPHIC  CF-GH 

The  development  of  Che  network  model  for  lexicographic  integer  GP 
probleas  follows  the  same  basic  procedure  as  outlined,  in  the  previous 
section,  for  weighted  integer  GP.  However,  we  shall  fora  a  sequence  of 
weighted  integer  GP-GN  models  for  the  representation  of  the  lexicographic 
model  where  the  total  number  of  models  in  each  sequence  shall  be,  at  worst, 
equal  to  the  number  of  priority  levels  in  the  original  lexicographic  GP 
mathematical  model. 

We  shall  denote  this  process  as  Sequential  GP-GN,  or  SGN.  SGN,  in 
turn,  is  simply  the  network  equivalent  of  the  approach  known  as  SGP  (l.e., 
sequential  goal  programming)  that  has  been  employed  for  approximately  15 
years  in  modeling  and  solving  large-scale  GP  models  [24,  27,  29].  The  pro¬ 
cedure  may  be  described  as  follows.  First,  we  formulate  that  portion  of  the 
lexicographic  GP  model  associated  with  priority  one.  This  submodel  is  simply 
a  weighted  integer  GP  model  and  thus  may  be  solved  directly  via  the  weighted 
GP-GN  approach  as  described  previously.  Next,  we  form  the  GP  model  associated 
with  priority  one  and  two.  However,  we  augment  this  model  with  the  additional 
condition  (i.e.,  a  new  rigid  constraint)  which  restricts  the  next  solution 
derived  to  one  which  does  not  degrade  the  solution  found  for  the  previous 
submodels  [i.e.,  g^(n,p)  *  a^*].  This  process  is  continued  until  all 
priority  levels  (and  associated  submodels)  have  been  solved.  The  last  solu¬ 
tion  is  the  solution  to  the  original  lexicographic  GP  model.  The  process 
may  be  summarized  via  the  algorithm  provided  below. 
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SCH  Algorithm 

Step  1:  Develop  the  lexicographic  integer  GP  model  as  shown  in  (2.5)  through 
(2.8).  Set  k  -  1. 

Step  2:  Establish  the  weighted  integer  GP  model  associated  with  priority 
k  (i.e.,  priority  one).  This  is: 
minimize  a^  «  g^  (n,p) 
s.t. 


ft(x)  +  ni  "  "  bi  1  £  P1 

n,  o  _>  0 

^  e  (0,1}  Vj 

Step  3:  Convert  the  present  GP  submodel  to  the  equivalent  weighted  integer 
GP-GN  model  and  solve,  via  the  appropriate  GN  software,  for  the 
optimal  value  of  a^  (i.e.,  a^*) . 

Step  4:  If  k  •  K,  stop  as  the  last  solution  (x*)  Is  the  solution  to  the 
original  lexicographic  GP  model.  If  k  <  K,  go  to  Step  5. 

Step  5:  Set  k  ■  k  +  1  and  establish  the  next  (weighted)  GP  submodel  in 
the  sequence  where  this  submodel  is  given  by: 
minimize  afc  -  gfe  (n,p) 
s.t. 


f.  (x)  + 


°i-bi 


i  e  T 


8t  (n,o)  -  at*  t-1,2, . . . ,k-l 

n,o  >  5 


e  {0,1}  VJ 

and  is  the  set  of  all  goals  and  constraints  in  priority  levels 


1  through  k 
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Step  6:  Repeat  Steps  3  through  5  until  the  antlra  nodal  has  baan  solvad. 

Tha  final  valua  of  x*  la  tha  optimal  program  for  tha  original 
lexicographic  integer  GP  modal. 

Implementation 

Implementation  of  tha  SGN  procaaa  would  be  achieved  via  the  use  of  a  buffer 
program  in  conjunction  with  the  GN  software.  The  buffer  program  simply 
automates  the  above  process  by  forming  each  submodel,  delivering  the  sub¬ 
model  to  the  GN  software,  and  forming  the  necessary  augmented  constraint 
associated  with  the  intermediate  GN  outputs.  In  actual  practice,  however, 
there  are  a  number  of  refinements  that  may  be  added  so  as  to  substantially 
improve  overall  computational  efficiency.  We  touch  on  just  a  few  of  these 
in  the  example  to  follow  as  well  as  in  the  next  section. 

Example 

In  order  to  simply  illustrate  the  basic  procedure,  we  shall  formulate 
the  sequence  of  models  associated  with  the  problem  shown  below: 

Find  x  to  lexicographically 
minimize: 

a  -  {(px),  (5n2  +  2p3),  (n3)}  (5.1) 


s.t 


xl  + 

x2  + 

X3  +  nl  -  °1 

-  2 

(5.2) 

3x 

1 

+ 

2x3  +  n2  "  °2 

■  4 

(5.3) 

3X1  + 

2x2  + 

7x3  +  n3  -  P3 

-  5 

(5.4) 

Xj  e  (0,1}  and  n,  p  >  0  (5.5) 

The  three  models  developed  to  support  the  SGN  process  are  shown  in  Figures 
5-1  through  5-3  and  we  shall  briefly  comment  on  the  construction  of  each. 
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In  Figure  5-1,  the  GP-GN  model  for  priority  level  on*  in  shown.  The 
solution  givss  a^*  *  0  ■  p  and  we  move  to  the  model  for  priority  levels 
one  and  two,  as  shown  in  Figure  5-2.  However,  we  may  simplify  constraint 
node  I  at  this  time.  That  is,  since  must  be  kept  to  a  value  of  zero, 
equation  (5.2)  is  equivalent  to: 

X1  +  x2  +  x3  +  nl  "  * 
or,  more  simply, 

X1  +  x2  +  x3  —  ^ 

As  a  result,  Che  goal  node  for  goal  I  (as  in  Figure  5-1)  may  be  replaced 
by  the  simpler  constraint  node  shown  in  Figure  5-2. 

Solving  the  GP-GN  model  of  Figure  5-2  gives: 

a2*  -  5  -  5n2  +  2p3  (5.6) 

which  is  the  augmented  constraint  that  should  be  added  to  the  next  GP-GN 
model.  However,  as  indicated  in  Figure  5-3,  relationship  (5.6)  can  be 
included  in  the  network  by  simply  modifying  the  two  goal  nodes  as  shown. 
That  is,  we  add  the  augmented  source  node  "A"  with  an  input  of  exactly  5 
units  (l.e.,  the  value  of  a2*).  The  two  arcs  emanating  from  node  A  then 
assure  us  that,  in  any  solution  to  the  GP-GN  network  of  Figure  5-3,  rela¬ 
tionship  (5.6)  must  be  satisfied. 

Solving  the  final  GP-GN  model  of  Figure  5-3  gives  us: 
a3*  «  0 

I  *  -  (1,1,0) 

and  thus  the  solution  to  the  example  of  (5.1)  through  (5.5)  is: 

I  *  -  (0,5,0) 
x  *  -  (1,1,0) 
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6.  REFINEMENTS:  LEXICOGRAPHIC  GP-GN 

Just  as  there  have  beea  found  to  be  nuaerous  refinements  to  SGP  (such 
as  "column  drop",  "early  stop",  etc.  [27]),  there  are  similar  refinements 
possible  for  SGN.  In  fact,  the  very  use  of  the  network  representation 
provides  for  a  number  of  straightforward,  but  potentially  powerful,  refine¬ 
ments  that  may  greatly  enhance  computational  efficiency.  We  discuss  a  few 
of  the  more  obvious  of  these  below. 


Augmented  Constraint 

Step  S  of  the  SGN  algorithm  indicates  the  addition,  for  each  submodel 
except  the  first,  of  the  augmented  constraint: 

gt  (n.p)  *  ac*  (6.1) 

From  this,  it  would  appear  that  we  are  increasing  the  complexity  of  the 
network  equivalent  when,  in  fact,  quite  often  just  the  opposite  occurs. 

To  illustrate,  consider  first  the  case  in  which  the  augmented  con¬ 
straint  (6.1)  contains  deviation  variables  from  more  than  one  goal  (i.e., 
there  exist  several  goals  at  the  previous  priority  level).  This  case  may 


be  considered  via  the  use  of  the  specific  example  given  below: 


gt(n,p)  -  +  u2n2 


(6.2) 


In  this  example,  we  have  two  goals  at  priority  t.  The  goal -nodes  for  these 
two  goals  are  given,  originally,  as  shown  in  Figure  6-1.  Note  that,  if 
the  achieved  level  of  (6.2)  is  zero  (i.e.,  a*  -  0) ,  a  simplification  is 
possible.  That  is,  if  a*  -  0  then  (6.2)  becomes: 


wlpl  +  u2n2  "  0 


which  can  only  be  satiofied  -./hen  there  is  no  flow  across  either  the  p^  or 
arc.  For  this  case,  we  simply  remove  these  two  arcs  from  their  goal-nodes 
of  Figure  6-1  in  all  further  calculations. 
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When  the  optimal  solution  to  (6.2)  Is  non-ssro  (i.e. ,  a*  >  0),  ws  ■ ay 
uss  tbs  nstvork  modification  slrssdy  illustrstsd  In  Figurs  5-3,  for  our 
example  of  ssetlon  5. 

Next,  1st  us  eonsidsr  tbs  csss  whsrsin,  for  ths  priority  level  previously 
solved ,  thers  is  but  ons  sssocistsd  gosl.  That  is,  ths  associated  augmented 
constraint  is  of  ths  general  fora  shown  below: 

gt(n,p)  -  wt  pt  +  u±  r>1  ■  a*  (6.3) 

Now,  if  the  value  of  a*  is  sero,  we  have  an  obvious  simplification  designated  as 
the  "collapse"  of  the  associated  goal  node.  This  is  illustrated,  for  a 
specific  fora  of  (6.3),  in  our  previous  example  of  section  5  wherein  the 
goal  node  of  Figure  5-1  "collapsed"  into  a  standard  constraint  node  in 
Figure  5-2. 

Variable  Reduction 

It  eay  also  be  possible,  when  performing  the  SGN  algorithm,  to  actually 
eliminate  one  or  more  variables  (and  thus  their  associated  nodes  and  arcs) 
from  the  model.  The  most  obvious  case  is  when  a  constraint  node  (either  an 
original  constraint  node  or  a  collapsed  goal-node)  can  only  be  satisfied  by 
one  specific  combination  of  those  variables  whose  arcs  lead  to  that  node. 

Early  Stop 

If,  through  variable  elimination  (as  directed  above),  all  the  variable 
values  are  fixed,  the  problem  is  solved  and  thus  it  is  not  necessary  to 
consider  further  submodels  in  the  sequence. 
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Honlinasr  Functions 

The  processes  previously  described  are  also  applicable  in  principle, 

to  nonlinear  integer  GP  models.  Here,  the  single  arc  from  a  variable  node 

to  a  constraint  or  goal-node  would  be  replaced  by  several  arcs.  Each  of 

these  arcs  would  represent  a  linear  approximation,  over  a  given  Interval, 

* 

of  the  nonlinear  function.  The  interval,  in  turn,  is  controlled  via  the 
upper  and  lower  bounds  on  flow  across  each  of  the  multiple  arcs  associated 
with  the  specific  variable-constraint/goal  combination. 

Mi mnax  GP  and  Fussy  Programing 

In  Section  2,  we  presented  only  two  (i.e.,  weighted  GP  and  lexico¬ 
graphic  GP)  forms  of  GP.  Another  form  of  GP  (see  the  references  for  further 
details  [27,  32,  41])  is  that  known  as  minmax  GP.  In  minmax  GP,  we  strive 
to  minimize  the  worst  single  goal  deviation.  Two  natural  and  straight¬ 
forward  extensions  of  minmax  GP  are:  lexicographic  minmax  GP  and  an  approach 
known  as  "fuzzy  programming"  [27,  32,  41].  By  means  of  approaches  similar 
to  those  discussed  in  this  paper,  these  models  may  also  be  represented  and 
solved  via  the  GN  approach  and,  consequently,  a  detailed  discussion  is 
believed  unnecessary. 

Augmented  GP 

Yet  a  further  class  of  GP  is  that  known  as  augmented  GP  [27,  34]. 
Augmented  GP  may  be  used  to  derive  (via  minimal  interaction  with  the  decision 
maker)  a  subset  of  all  the  nondomlnated  solutions  which  exist  for  a  given 
problem.  Such  an  approach  is  Implemented  by  means  of  extensions  to  the  GP-GN 
concept  and  is  the  topic  of  a  present  research  effort. 

*The  restrictions  typical  of  all  such  approximations  must,  of  course,  hold. 


7.  STOMUH 


In  this  paper  v*  have  providad  and  dascrlbad ,  via  example,  tha  raaulCs 
of  an  on-going  raaaarch  afforc  dadicatad  to  tha  nodal log,  solution  and 
analysis  of  nultiobjactlva  intagar  prograanlng  problana.  This  has  baan 
acconplishad  by  naans  of  combining  ganeralizad  goal  programing  [21,  24, 

25,  27,  28,  29,  30,  34}  with  generalized  networks  11,  2,  3,  4,  5,  6,  7,  8, 
9,  10,  11,  12,  13,  14,  15,  16,  17,  18,  19,  22,  26,  31,  33,  35,  37,  38,  39, 
40].  Tha  resulting  methodology,  denoted  as  GP-GN  may  be  applied  to  any 
multiple  objective  model  which  conforms  to  the  standard  assumptions  of 
generalized  networks.  For  such  problems,  the  GP-GN  approach  should  provide 
for  a  computationally  efficient  method  of  problem  solving. 
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